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We extend our previous analysis on the exchange statistics of vortices having a single Dirac fermion 
trapped in each core, to the case where vortices trap two Dirac fermions with U(2) symmetry. Such 
a system of vortices with non-Abelian Dirac fermions appears in color superconductors at extremely 
high densities, and in supersymmetric QCD. We show that the exchange of two vortices having 
doublet Dirac fermions in each core is expressed by non-Abelian representations of a braid group, 
which is explicitly verified in the matrix representation of the exchange operators when the number 
of vortices is up to four. We find that the result contains the matrices previously obtained for 
the vortices with a single Dirac fermion in each core as a special case. The whole braid group 
does not immediately imply non-Abelian statistics of identical particles because it also contains 
exchanges between vortices with different numbers of Dirac fermions. However, we find that it does 
contain, as its subgroup, a genuine non-Abelian statistics for the exchange of the identical particles, 
that is, vortices with the same number of Dirac fermions. This result is surprising compared with 
conventional understanding because all Dirac fermions are defined locally at each vortex, unlike the 
case of Majorana fermions for which Dirac fermions are defined non-locally by Majorana fermions 
located at two spatially separated vortices. 

PACS numbers: 05.30.Pr, 74.25. Uv, 67.85.-d, 21.65. Qr, 03.67.-a 



I. INTRODUCTION 

Topological insulators/superconductors have an attractive property that some of them possess quantum vortices 
which trap zero-energy, Majorana or Dirac, fermions in their cores m,l3- The existence of such zero-energy fermions is 
topologically protected and is robust against small perturbations Thus, when we consider adiabatic manipulation 

of vortices such as interchanging the positions of two vortices, we can treat the vortices as objects that are always 
accompanied by zero-energy fermions. In particular, according to the recent discoveries, the exchange of such vortices 
can be represented by a non-trivial representation of a braid group, whose precise form is determined by the trapped 
zero-energy Majorana fermions [sl-fioj and Dirac fermions [ll| . The exchange of vortices with Majorana fermions gives 
a non-Abelian statistics because they are all identical particles. The statistics is called non-Abelian because quantum 
states of two vortices transform non-diagonally under the exchange of two vortices (i.e., the exchange operation is 
described by non-diagonal matrices acting on the quantum states), and two adjacent exchange operations (such as 
those for vortex pairs (1,2) and (2,3)) do not commute with each other. This is highly contrasted with the ordinary 
statistics where only a phase factor e*^ appears under the exchange of two particles {9 = for the Bose-Einstein, 
6 ~ TT for the Fermi-Dirac, and others for the anyon statistics), and two adjacent operations arc commutative. On the 
other hand, it is unclear whether the exchange of vortices with Dirac fermions gives a non-Abelian statistics or not 
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because, in general, it exchanges different particles, that is, vortices with different numbers of Dirac fcrmions. For 
U(l) Dirac fermions, vortices can be distinguished by the occupancy of Dirac fermions. 

According to the classification of topological insulators/superconductors [H, 0] and its extension to the case with 
topological defects the vortices with the Majorana or Dirac fermions are categorized into different types: class 

D for the Majorana and classes C and Dili for the Dirac. However, the essential difference between the vortices with 
Majorana fermions and Dirac fermions is the parity of the number of zero-energy Majorana fermions trapped to a 
single vortex. Notice that a single Dirac fcrmion corresponds to two Majorana fermions. Thus, when a vortex traps an 
even number of the Majorana fermions, it should be regarded as the Dirac case, while an odd number of the Majorana 
fermions, just as the Majorana case [1]. So far, non-Abelian unitary transformations (non-Abelian representations of 
a braid group) are found both in the Majorana and Dirac cases. In the Majorana case, the non-Abelian statistics was 
first discovered when a vortex has a single Majorana fermion @ and later when a vortex can have multiple Majorana 
fermions with non-Abelian symmetry j8l4lOl|. In contrast, in the Dirac case, while non-Abelian representation of a 
braid group is analyzed when a vortex traps a single Dirac fermion [llj , it is unclear if it gives a genuine non-Abelian 
statistics. [23 The present paper discusses the non-Abelian representation of a braid group for the exchange of vortices 
which have multiple Dirac fermions with non-Abelian symmetry, to complete the series of analyses. We then show 
that it indeed contains a genuine non-Abelian statistics of the exchanges of identical particles, that is, vortices with 
the same numbers of Dirac fermions. This result is somewhat surprising compared with conventional understanding 
because non-Abelian statistics appears in a system with only locally-defined Dirac fcrmions. 

There is an important difference between the vortices with Majorana fermions and Dirac fermions. Consider a 
system of vortices each of which traps only a single Majorana fermion. Recall that the Majorana fermion has a 
unique property that there is no distinction between a particle and a hole (anti-particle) [Tsj . Thus, in order to define 
the Fock vacuum, one has to introduce a Dirac fermion by using two Majorana fermions that belong to different 
vortices [1, Q. Therefore, the Dirac fermions thus constructed are non-local objects. When there are 2m vortices, 
the total Hilbert space of the zero-energy fermions has a dimension (degeneracy) of 2™, where each zero-energy Dirac 
fermion has the dimension 2 (empty or occupied) . On the other hand, in the system of vortices each of which traps a 
sing le Dirac fermion, we can immediately construct the Hilbert space without introducing 'non-local' Dirac fermions 
, and find a non-Abelian representation of a braid group (but not non-Abelian statistics) for the exchange of two 
vortices. When there are 2m vortices, the dimension of the Hilbert space is 2^™ = 4™. 

In the Majorana case, increasing the number of Majorana fermions in a single vortex brings in an interesting 
nontrivial structure. When the multiple Majorana fermions are in the vector representations of SO (3) @ and, in 
general, S0(2iV - 1) the exchange matrices of two vortices are given as tensor products of the matrices that 
appear in the single Majorana case (called the Ivanov matrices) and generators of the Coxeter group of the ^2m-i 
type (for 2m vortices). Besides, the dimension of the Hilbert space of zero-energy fermions becomes larger than that 
of the single Majorana case by the internal degrees of freedom. For example, when three (triplet) Majorana fermions 
with SO (3) symmetry are trapped, the dimension is 2"^™ which should be compared with 2™ for the case with a 
single Majorana fermion. It is not known if a similar structure exists and how the dimension of the Hilbert space is 
enlarged, when increasing the number of Dirac fermions at each vortex. The purpose of the present paper is to show 
the explicit forms of the non-Abelian representation of a braid group for vortices with multiple Dirac fermions, and 
that it contains a genuine non-Abelian statistics in the sectors of the exchange of vortices with the same numbers 
of Dirac fcrmions. As the simplest but non-trivial example, we focus on the Dirac fcrmions with U(2) symmetry. 
Extension to general cases will become more complicated, but should be straightforward. 

It should be noticed that, since the vortices with zero-energy fermions are characterized by topology, they appear 
in many different quantum systems. For example, vortices with the Majorana fermions are realized in chiral p-wetve 
superconductors Q , such as Sr2Ru04 , in chiral p-wave superfiuids, such as the A-phase of '^He in 2-|-l dimensions, 
and also in other systems [l6j. Vortices with the Dirac fermions are identified with the integer (singular) vortices 
in the '^He A-phase in 2-1-1 dimensions 17| the normal "o" vortices in the '^He B-phase in 3-1-1 dimensions [l^ and 
also with dislocation lines in topological insulators [l^. In these examples, the vortex has only a single zero-energy 
Majorana or Dirac fermion at its core. So far, there is no condensed-matter example of vortices with more than 
two Majorana fermions. However, we know at least two examples in high-energy physics. In fact, the existence of 
zero-energy modes in the vortex-fermion system was discussed long time ago in the context of relativistic quantum 
field theory 0. 

The primary example is the color superconductor in QCD which could exist in extremely high density matter such 
as in the cores of neutron stars [2l| . In particular, vortices with non-Abelian symmetries appear in the color-flavor 
locked (CFL) phase where the original color SU(3) and flavor SU(3) symmetry in the vacuum breaks down to SU(3) 
with color and flavor degrees are locked (2^ . The SU(3) CFL symmetry is further broken down to a U(2) symmetry 
in the core of non-Abelian vortices [2^ . It has been recently shown [23| by using the index theorem that there appear 
two types of non-Abelian vortices: the vortex which traps a triplet of zero-energy Majorana fermions Q and the 
vortex which traps an U(2) doublet of zero-energy Dirac fermions. While the former gives an example of non-Abelian 
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statistics of vortices with multiple Majorana fcrmions the latter gives an example of non-Abelian representation 
of a braid group for the exchange of vortices with doublet Dirac fermions focused in the present pap er. The secondary 
example is the non-Abelian vortices in supersymmetric QCD with U(A^) gauge symmetry [25| (see Ref. [1^ for 
reviews). In this case, the color U(A'^) symmetry and the flavor SU(-/V) symmetry are spontaneously broken down to 
the SU(Af) color-flavor locked symmetry by the scalar quark condensates in the vacuum, and it is further broken down 
to a U(iV — 1) symmetry in the core of non-Abelian vortices. These vortices contain one singlet and one {N — l)-plet 
of zero-energy Dirac fcrmions in their core. Therefore, it gives an example of non-Abelian statistics of vortices with 
arbitrary number of Dirac fermions. 

Before finishing Introduction, let us briefly comment on the potential application to quantum computers [13, . 
As discussed above, the system of vortices with zero-energy fermions is robust against small perturbations from 
environment, and has the Hilbert space with a large dimension. These are desirable properties as quantum computers. 
Comparing the Majorana and Dirac cases, it should be noticed that the Dirac case is simpler because we do not have 
to introduce non-local Dirac fermions. Since vortex systems with non-Abelian symmetry have larger dimensions, 
it would be worth to consider the case with non-Abelian Dirac fermions, even though it has not been realized in 
laboratory. 

This paper is organized as follows. In Section |lll we briefly summarize the non-Abelian statistics for the vortices 
trapping a single Dirac fermion with U(l) symmetry at each core, as presented in Ref. [llj. In Section Hill wc discuss 
the non-Abelian representation of the braid group for the exchange of the vortices trapping doublet Dirac fermions 
with U(2) symmetry. In Section HVl we discuss the difference between non-Abelian representations of the braid group 
for the exchange of the Dirac fermion with U(l) symmetry and that with U(2) symmetry. In Section |Vl we show 
that the whole braid group contains a genuine non-Abelian statistics as its subgroup, in the sectors of the vortices 
with the same number of Dirac fermions. Section IVII is devoted to a summary. In Appendices, we present detailed 
supplementary information. In Appendix |X1 we give the transformation matrices for U(l) (singlet) Dirac fermions. 
In Appendix IB] we give the Hilbert space and the exchange matrices for n = 4 U(2) Dirac vortices. In Appendix [Cl 
we discuss restricted Hilbert subspaces, in which two successive exchanges of vortices is equal to the identity. 



II. NON-ABELIAN REPRESENTATION OF BRAID GROUP FOR EXCHANGE OF U(l) DIRAC 

VORTICES 

Let us first explain how the non-Abelian representation of the braid group appears in the system of vortices having 
a single Dirac fermion in each core, which corresponds to considering the Dirac fermion with the U(l) symmetry. This 
is a brief summary of the recent work done by three of us [ll| . We highlight the similarities to and differences from 
the case with Majorana fermions. Below, we call vortices with the Dirac fermions "the Dirac vortices", while vortices 
with the Majorana fermions, "the Majorana vortices" . 

Consider n Dirac vortices which are labelled by k ~ 1, • • • ,n. The number of vortices, n, can be arbitrary in 
contrast with the case of Majorana vortices where we define Dirac fermions by using two Majorana fermions, and 
thus the total fermion number is even. The operator denotes the Dirac fermion of the A:-th vortex. Together with 
its hermitian conjugate tpj,, they satisfy the following algebra: 

{^k,i^l} = Ski, {i^k,M = o. Wi,^l} = o. (2.1) 

We regard tj^k and tpj. as the annihilation and creation operators, respectively. Exchange of the k-th and (fc + l)-th 
vortices, which is denoted by T^, induces the following exchange of ipk and V'fe+i, 

T,:\f ^^f' , (2.2) 

L i'k+i -Vk 

with the rest i/if (^ 7^ and /c + 1) imchanged. We note that T^'s satisfy the braid relations, 

(i) TkTiTk ^ TiTkTe for \k - £\ = 1 , (2.3) 

(ii) TkTi = Ten for \k^e\>l, (2.4) 

as a general rule of exchange operations. Notice that these relations and the transformation (|2.2p are the same as in 
the case of the Majorana vortices, while the fermion operators tpk and '0^ now satisfy Eq. (|2.ip . instead of the Clifford 
algebra {7fc, ji} = 2Ske satisfied by the Majorana fermion operators 7fc's. Therefore, one has the same property of the 
operator Tk, namely, four-time successive application of is equivalent to the identity, (T^)* = 1. This fact by itself 
suggests that the exchange of these vortices shows a representation of the braid group different from the ordinary 
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Bose-Einstein, Fcrmi-Dirac or Abelian anyon one which shows Abehan representation of the braid group. In fact, 
the exchange statistics of the Majorana vortices showing the same property (T^)^ = 1 turned out to be non- Abehan 
0, and we will discuss below if the same is true for the Dirac vortices. 

The transformation (|2.2p can be represented by the following unitary operator 

tI = i + ^fc+iV-l + ^l+,i^k - V-i+iV^fc+i - i^li^k + s^Al+iV^fc+iV^l^- , (2.5) 

so that T^'ipi{T^)^^ {£ ~ I,-- - ,n) reproduces the transformation law. The superscript "s" implies the singlet Dirac 
fermion. One can also confirm by a straightforward calculation that f^'s satisfy the braid relations, 

(i') f^-f|f| = f|f|f| for |fc-£| = l, (2.6) 
(ii') f|f| = f|f| for |fc-^|>l. (2.7) 

Having the explicit form of the exchange operator (|2.5p . we are able to check the representation of the braid group 
for the exchange of the Dirac vortices. First of all, as we mentioned above, four successive exchanges of the fc-th and 
(fc + l)-th vortices indeed yield the identity: 

inr = h (2.8) 

while two successive exchanges do not, 

{fiy = (1 - 2^^4)(1 - 2^l_^,^,^,) + 1 . (2.9) 

Next, it should be noted that f| and t\j^^ are not commutative; 

[T|,r|+i]^0. (2.10) 

Therefore, one observes at the operator level that the exchange operation of two Dirac vortices is in general non- 
Abelian. However, this observation does not necessarily mean that the exchange is always non-Abelian. In order to 
confirm whether the exchange operation is indeed non-Abelian, we have to check if the commutator [f^ , t^j^-^ does 
not vanish in the matrix representation. 

Let us consider the matrix representation of the operator f| . As a basis of the Hilbert space, we choose the Fock 

states defined by the Dirac fermion operators V'fc's. One of the merits of the Dirac vortices is that we can construct 
the Hilbert space naturally by using locally- defined Dirac fermions, which is in clear contrast with the Majorana case. 
We first define the Fock vacuum state |0) by 

V)f|0) = for aU i. (2.11) 

Then, by acting successively the creation operators -i^j's on the vacuum, we obtain the /-particle state {0 < f < n) 

I0...01...1...10---0) =Vi,\...7A,V..ViJjO), (2.12) 
i ii h if n 

in which the ^^-th (i = 1, • • • , /) vortex is occupied by a Dirac fermion, while the other vortices are empty. 

When we have only one vortex, n = 1 (fc-th vortex), there are two Fock states, |0) and |1) = ■(/'|.|0). 

When we have two vortices, n = 2 (fc-th and (fc + l)-th vortices), there are 2^ = 4 Fock states, |00) = |0), 
|10) ^ V^IlO), |01) ^ V^I+JO) and |11) ^ ^I^I+JO). 

We can similarly obtain the Fock states for any number of vortices n. The basis of the whole Hilbert space for the 
n- vortex system is given by a tensor product of the 4 Fock states constructed at each vortex. Because the fermion 
number operator P = '^"^i'4'l'4'i commutes with f| for £ = 1, ■ ■ ■ ,n — 1, the fermion number / (an eigenvalue 
of /^) is a conserved quantity under the exchange of vortices. Thus, the whole Hilbert space H^"^ for n vortices 
(with the dimension 2") can be decomposed into subspaces M^"'^^ labelled further by the total fermion number /: 
H(") = e^^-^H^"^-''). Then, in each subspace, the operators f|'s are represented by matrices, whose explicit expressions 
are shown up to n ~ 4 in Appendix [X] One can confirm that the matrix representations of the exchange operators 
ff's are indeed non-Abelian for n > 3. 
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III. NON-ABELIAN REPRESENTATION OF BRAID GROUP FOR EXCHANGE OF U(2) DIRAC 

VORTICES 

Now let us turn to the case with vortices which trap two massless Dirac fermions having the "pseudo-spin" U(2) 
symmetry. In particular, we consider Dirac fermions in the doublet of U(2), which are denoted by (« = I7 2) for 

the A:-th vortex. We also use the vector notation for the doublet {tpl,tplY. The exchange of the fc-th and (fc + l)-th 
vortices induces transformation of the Dirac fermions. Here we consider the simplest transformation similarly as the 
single Dirac fcrmion shown in Eq. (j2.2p : 



-* -*% 

for each a = 1,2, with the rest ipf {£ k and fc + 1) unchanged. This transformation is expressed by the unitary 
operator 

^ n (1 + + ^tiir^ - ^tik+i - r^Hi + 2V',ti^fe+iVif vifc) , (3.2) 

a=l,2 

which is invariant under the U(2) transformation, (V'feiV'fe)* ~^ exTp{iip) exTp{id ■ (t/2)(?/;^,, -0^)* with ip and 6 being 
parameters and a the Pauli matrices. We confirm that ff^ijjgT^^ (a = 1, 2) reproduces the transformation p.ip . and 
that ffc's satisfy the braid relations, (i') ffc-r^ffc = fifkTe for |fc — £| = 1 and (ii') Tkfi = titu for |fc — £| > 1. Similarly 
as in the U(l) case (see Eqs. (|2.8p and p.9p ). the exchange operator satisfies 

= (3.3) 

while two successive exchanges do not go back to the identity at the operator level: 

(ffc)2 = (1 - 2^im){l - 2^mi){l - 2^ll,^l+,){l - 2^ll^i^l_,^) ^ 1 . (3.4) 

We can also check that ffe and Tk+i are non-commutative, which suggests that the representation of the braid group 
for the exchange of U(2) Dirac vortices is non-Abelian. 
We define the number operator of the Dirac fermions 

n 

£=1 a=l,2 

whose eigenvalues / give the total number of zero-energy Dirac fermions in the vortex system (0 < / < 2n). We 
note that X]a=i ^''Pk^^k U(2) invariant, hence / is also U(2) invariant. Wc also note that / commutes with fg for 
^=1,... ,n— 1, hence / is a conserved quantity under the exchange of two Dirac vortices. 

We can construct the Hilbert space by successively applying -0°^ to the Fock vacuum |0) defined by 

V)"|0) = for aU I and a = 1, 2. (3.6) 

Below, we explicitly construct the Hilbert spaces when the numbers of vortices are n = 1, 2, 3, and see the non-Abelian 
properties of the exchange operation in the matrix representations for n = 2 and 3. The results for n = 4 vortices are 
quite involved and are shown in Appendix IB] 



A. The case of n = 1 



When we have only one vortex, we are not able to discuss the exchange of vortices. However, let us consider this 
case to demonstrate how to construct the Hilbert space for the U(2) Dirac vortex. We have two massless Dirac fermion 
operators ij}], and ifr^ (to avoid notational confusion, we use the label k to specify this single vortex). By applying the 
creation operators -ip]^ and to the Fock vacuum |0), we obtain 2x2 = 4 energetically degenerate states (2 for 
empty/occupied, and another 2 for a = 1,2). Let us introduce the notation \'R-f) where TZ is the representation of 
U(2) group and / is the total number of the Dirac fermions. Namely, 



|lo) ^ |0>, 
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|2i) ^ (|I)|0)' (3-7) 

The bold numbers 1, 2 imply the singlet and doublet representations of the U(2) group. Notice that the fully occupied 
state ^fc^^fc^|0) is invariant under the U(2) transformation, thus it belongs to the singlet representation. Hence, there 
are two singlet states |lo) (empty) and II2) (fully-occupied by two fcrmions), and one doublet state |2i) (occupied by 
one fermion). We use these states as the basis to span the whole Hilbert space of the zero-energy states. Therefore, 
we have decomposed the representations of the U(2) pseudo-spin as 

I0 + 21 + I2, (3.8) 

and correspondingly the whole Hilbert space of a single vortex H^"=^^ into a direct sum 

= H^m e H^w ® M^i^i , (3.9) 

where H^m = {|lo)}, H^w = {|2i)} and H^i^i = {|l2>}. 



B. The case of n = 2 



Consider the case when we have two Dirae vortices which are respectively labelled by k and fc -I- 1. First of all, we 
can use the decomposition for a single Dirac vortex, Eq. p.Sp . to decompose the whole Hilbert space of two Dirac 
vortices into representations. Since each vortex contains the representations (Iq + 2i + l2)vortexf = k, k + 1), the 
total representation is obtained as a tensor product of them, which can be decomposed as 

(Iq + 2i + l2)vortcxfc ^ (Iq + 2i + l2)vortexfe+l 

= loo + 111 + I20 + I02 + I22 

+ 2io + 2oi + 221 + 2i2 

+ 3n, (3.10) 

where the bold numbers denote representations, and the subscript UkUk+i {uk, nk+i = 0, 1, 2) denotes the number 
of the Dirac fermions, Uk and rife+i, at the fc-th and {k + l)-th vortices, respectively. Next, one obtains the basis of 
the whole Hilbert space by applying ■0"^ {£ ~ k, k + 1 and = 1,2) successively to the Fock vacuum |0) defined 
by -0^1 0) =0 for all i and a = 1, 2. Then, one reorganizes the Fock states according to the decomposition into 
representations . 

Below, we explicitly show the basis of the Hilbert space according to the decomposition in Eq. (|3.10p . We introduce 
the notations iTZn^nk+i) for basis states and H'^i^i^ai for the subspaces of the Hilbert space. Here TZ denotes the 
representation of pseudo-spin, and the subscript UkUk+i {nk, rik+i = 0, 1, 2) denotes the number of the Dirac fermions 
at the fc-th and (fc -I- l)-th vortices, respectively. In the subspaces, we do not distinguish Uk and Uk+i for the reasons 
mentioned later, and thus iVi and iV2 are defined as A^i = max{nfe,nfc+i} and N2 = min{nk,nk+i}- For example, 
(TVi, A'2) = (1,0) contains two cases (rifc, rifc+i) = (1, 0), (0, 1). Then, the whole Hilbert space of two Dirac vortices is 
decomposed into seven subspaces: 

Jj{,l=2} ^ Jjl[00J ^ Jjl[llJ ^ Jjlpo] ]e[l[22] 

® H^["i . (3.11) 

This should be compared with the decomposition (|3.10p . 

The first line in Eq. (|3.1ip corresponds to the singlet subspaces. They are defined by the basis states as M^i°°'i = 
{|loo)}, H^im = {|lii)}, H^poi = {|l2o), II02)}, and H^i^^i = {|l22)}- Explicit forms of the basis states are given as 
follows: 



|loo) ^ |0), 
|lii> 

|l20> 



^(^r^fcii-^r^iii)io), 

ltj,2t 



^k'rk'm, 



(3.12) 
(3.13) 

(3.14) 
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|lo2) ^ ^llii^lU^), (3.15) 
II22) ^ ^I^V'fV'^i.i^fc^ilO). (3.16) 

The second line in Eq. p. lip corresponds to the doublet subspaces. They are defined by the basis states as 
p[ioi = {|2io), |2oi)} and H^pii = {|22i), |2i2)}. Explicit forms of the basis states are given as follows: 

|2io) ^ (fi)\0), 



k 



'Pk+l 
'Pk+l 



|0), (3.17) 



1221) ^ jio>, 

'^k \ „?,lt J,2t 



V^Vi^'fcVilO). (3.18) 



Here, in the last two expressions, we have factored out the singlet part so that the vector structure becomes manifest. 

Lastly, the third line in Eq. (|3.1ip corresponds to the triplet subspace. It is defined by the basis states as 
{|3ii)} with 



3[lll — 



>ii/ = I 72^^k^i\i + '^ki^lli) I |0)- (3.19) 

Now we have prepared to discuss the exchange of vortices. We recall that the exchange of ■0^ and "(AjJ^^ expressed 
by the operator ffc defined in Eq. p.2p as a unitary transformation Tk'il'iT^^ (£ = fc, fc + 1 and a = 1, 2). With the 
explicit forms of the basis states in Hilbcrt subspaces, we are able to express the operator fk as matrices. In the 
singlet subspaces H-*-!""!, H-*-!"!, M^i™^ and H-'-i^^i, the corresponding matrices are 

^Moo] _ -, 
^k " -'^i 

^k ^ 

1[20] _ f 1 

'^U — 



1 



r-- = 1. 

In the doublet subspaces M^[i"i and H'^i^ii, the corresponding matrices are 

2[io] _ / -1 



~ U 

i[21] _ / 1 

^fe " I -1 



Lastly, in the triplet subspace H^ini, the matrix is 



3r: 



Since H-'-tooi, H-*-!"!, H-'-P^i in the singlet subspace and H^i"! in the triplet subspace consist of only one basis state, we 
have one-dimensional representations. In contrast, the other subspaces H-'-Poi, H^ii"!, and H'^pii have two basis states, 
thus yielding two-dimensional representations. Notice that the exchange matrices in these subspaces have off-diagonal 
elements. However, this does not mean that the representation of the braid group is non-Abelian. Rather, this simply 
implies that two basis states are mixed with each other by the exchange operation (this is the reason why we did not 
specify the order of Uk and Uk+i in defining the subspaces). In fact, one can choose appropriate basis states so that 
the exchange matrices are expressed as diagonal. For example, Tj,'^"' is diagonalized with eigenvalues ±1 showing the 

Abelian representation of the braid group, while t^'^"' and t^'^^' arc diagonalized with ±i, showing anyon-like Abelian 
representation of the braid group. 



8 



C. The case of n = 3: emergence of non-Abelian representation 



Let us finally consider the case of three Dirac vortices which are respectively labelled by fc, k + l and fc + 2. Basically 
we will follow the same procedures presented in the case of two Dirac vortices, but as easily expected, the analysis 
becomes quite involved. Still, we present here all the information since this is the simplest case where the non-Abelian 
representation of the braid group appears. In fact, when we have only two vortices (fc-th and (fc + l)-th vortices), there 
is only one exchange operator f^., thus we are not able to discuss the non-commutativity of two exchange operations. 
It makes sense only when we have three or more vortices. When we have three vortices (fc-th, (fc + l)-th and {k + 2)-th 
vortices), we can define two exchange operators ffe and fk+i, and thus we can check if these two are commutative or 
not. 

First of all, by using Eq. (|3.8p for a single Dirac vortex, we can decompose the pseudo-spin structure made of three 



similar to the 



2 case. 



Dirac vortices into representations of U(2). We again introduce the notation TZ 

Here the subscript nknk+ink+2 [nk, '^-fc+i, = 0, 1, 2) denotes the number of Dirac fcrmions at the fc-th, (fc-l-l)-th 
and (fc + 2)-th vortices, respectively. Then, one finds 



(lo + 2i + l2)vortcxA; ® (lo + 2i + I2 ) vortex fc+1 ® (lo + 2i + l2)vortcxfc+2 
= loOO + (I02O + I2OO + I002) + (lllO + loil ^ 
+ (2l00 + 2oiO + 2ool) + (2210 + 2i20 + 2oi2 ^ 

+ (3iio + 3oii + 3ioi) + (3ii2 + 8211 + 3i2i) 
+ 4s 111 . 



lioi) + (I220 + I022 + I202) + (I112 + I211 + I121) H 
2io2 + 2o2i + 2201) + (2212 + 2i22 + 2221) + (2aiii 



1222 
" 2s 111) 



(3.20) 



Here we have also introduced new notations 2aiii, 2s m and 4s m, meaning that the first two subscripts with 
underlines are made Asymmetric (Symmetric) with respect to the indices. We will note again when we present the 
explicit forms of the basis states. Each term corresponds to the basis state in the Hilbert subspace. Notice that we 
have already grouped the representations so as not to distinguish the ordering of n^. rik+i and nk+2- Thus, one can 
easily see that the whole Hilbert space H{"=3} of three U(2) Dirac vortices can be decomposed into a direct sum 



Lpil] 



'-[222] 



j2[ioo] ; 

i3[iio] : 

j4[iii] 



[[^[2101 ] 
[T3[211] 



'[111] 



/at 



(3.21) 

fc, fc+1, fc-|-2 and a = 1, 2) successively 



Next, one can obtain the basis states of the Hilbert space by applying ij:^ 

to the vacuum |0) defined by V'^ |0) =0 for all i and a = 1, 2. 

The first line in the decomposition p.21[) corresponds to the singlet subspaccs. They are defined by the basis states as 

Hlln™l = {|looo)}, tfl^nnl = { | lo20> , I W , 1 1002) } , H^I""! = { | Ino) , | Iqii) , | llOl) }, ^P™! = { 1 1220> , 1 1022) , 1 1202) }, 



H^piii = {II112), II211), |li2i)} and H^P22] ^ {|l222)}- Notice that ah these states have even numbers of Dirac 
fermions. Explicit forms of the basis states are given as follows: 



l-[200] 



MllO] 



>-[220] 



'-[211] 



^V-i^iV^f+ilO), 

= V^fe+iV^fe+iV^fc+zV'fcUlo)' 



-•-000/ 

1020) 
1200) 
1002) 

liio) 
Ion) 
lioi) 

I220) 
I022) 
I202) 

I112) 



)|0), 



i^mk+i)i^l+2i'k+2i 



71 '/'fe^V'^ (V'fe+lV'fe+2 - V'fe+l'/'fc+2 



(3.22) 
(3.23) 

(3.24) 

(3.25) 

(3.26) 
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•■[222] 



/2t 



It ,72t 



(3.27) 



It is interesting to notice that the structures in H^i""! and H^i^ni, or in H^POQi and H^P^oi are similar. This is the 
symmetry between "occupied" and "empty" states (or "particles" and "holes"), reflecting the ambiguity in defining 
the creation operator (i.e., we could define tp^ as a creation operator, instead of the annihilation operator). 

The second line of Eq. p.2ip corresponds to the doublet subspaces. These four subspaces are respectively defined 
by the basis states as H^dooj = {|2ioo), |2oio), |2ooi)}, H^Pioi = {|22io), |2i2o), |2oi2>, |2io2>, |2o2i>, |22oi)}, IH[2[22i] = 
{I2212), |2i22), |222i)}, and H^ii"! = {|2aiii), |2sjLii)}. This time, the total numbers of Dirac fermions are odd in 
these states. Explicit forms of the basis states are given as follows: 



^[100] 



^[210] 



^[221] 



^[111] 



2ioo) — 



'010; 



2ooi) — 



-120) 



'012/ 



'021/ 



'201/ 



2122) — 



2a 



111, 




'S 111, 



2t 
+2 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



|0) 



Now the meaning of the notations A and S is evident. For example, in |2a 111) the first (fc-th) and second ((fc + l)-th) 
vortices form anti-symmetric combination with respect to the indices. 

The third line of the decomposition p.2ip corresponds to the triplet subspaces. They are defined by the basis states 
as H^'i""! = {|3iio), |3oii), |3ioi)} and H^i^ni = {|3ii2), |32ii), |3i2i)}. Total fermion numbers are now even. Explicit 
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forms of the basis states are given as follows: 



i[110] 



i3no>^ I ^M^V+i + r,'r^+i) I io) 



>oii; 



|3ioi) 



it 



+ ttlitlU) I |o), 



(3.32) 



fk+ltk+2 

timu 



Vk rk+2 



5[211] 



V2 



|321l) ^ tlHl^ 

I 



( 



Vk Vk+1 



|0), 



It 



rk+lVk+2 

-2t 



'It 



V 



|0), 



(3.33) 



V_k Vk+2 



Vk+lVk+2 



WkHV+2 + i^V^]}+2) UfcliV-^lilo). 



V 



72 



Vk 'Pk+2 



The last line of the decomposition p.2ip corresponds to the quartet subspace. It is defined by the basis states as 
I^dii] = {|4siii)}. Explicitly, 



i[iii] 



|4siii) 



f;Vk Vk+l^k+2 ^ J^\^k Vk+1 ^ Yk Vk+1 



V3 



^K'Vk Yk+l 

7=M^tlli + tVtlUtlU + ^tl^tlWkU 



)tll2 



V 



|0). 



Now we can finally discuss the exchange of U(2) Dirac vortices in the Hilbert space constructed for three Dirac 
vortices. As in the case of n = 2, we express the operators fj. and Tk+i as matrices with the basis states presented 
above. 

H[i[2oo]^ H^[ii«i, H^P^oi, M^[2iii and H^i^^^i, the exchange matrices are 

•■[000] 



In the singlet subspaces H^io""!, 



'-[200] 



L[1101 



■■[220] 



L[211] 



•■[222] 



[000] 

+1 


= 1, 




1 






1 


;) 











-1 














1 










1 









!) 


. 


1 






-1 














1 





-1 






J-[200] 
'fc+1 



-l[110] 

'fc+1 



J- [220] _ 



-l[211] 
'fc+1 




1[222] 
'fc+1 



(3.34) 



In the doublet subspaces, 



2[ioo]^ H^i^if'i, ]EII^p2i] and H^imi, the exchange matrices are 



'[100] 




^[100] 

'fc+1 
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/O 
1 



2[210] 



0-1 




2[iii] 



In the triplet subspaces, H^i""! and M^i^"!, the exchange matrices are 



[110] 



i[211] 




■^[110] 
'fc + 1 



■i[211] 
'fe + 1 



Lastly, in the quartet subspace M^imi, the exchange matrices are 



i[iii] 



«[iii] 
'fe+i 



1. 




1\ 



(3.35) 



(3.36) 



(3.37) 



We now find that, except for the onc-dimcnsional representations (in IHI-'-i'""'! , H-'-P^^i and IHI'*[i"i), all the exchange 
matrices are non-commutative; [Tk,Tk+i] 7^ 0. Therefore, the exchange operations in these subspaces exhibit non- 
Abelian representation of the braid group. 

Recall that we have shown that the exchange operators f^'s satisfy (f^)^ = 1 at the operator level (see Eq. p.3p ). 
Of course this is satisfied by all the exchange matrices obtained above. However, interestingly, a stronger relation 



[220] and H^iiiiJ satisfy this stronger relation. More 



(f^)^ = 1 is satisfied in some subspaces. In addition to the trivial one-dimensional representations in H-"-!"'"'! , H-'-P^^i 
and ]HI'*[iiii, we observe that the exchange matrices in H-'-i^ooi , 
details are discussed in Appendix [Cl 

We can continue the construction of the Hilbcrt space when the number of vortices is more than three. The 
dimension of the total Hilbert space of n-vortex system is (2^)". We present the results for n = 4 vortices in 
Appendix [BJ which also show non-Abelian representation of the braid group. 



IV. U(l) DIRAC STRUCTURE EMBEDDED IN U(2) DIRAC VORTICES 

Recall that we already found in the previous paper [ll| that the system of U(l) Dirac vortices shows non-Abelian 
representation of the braid group. Then, it is natural to ask how the previous results are related to the present results. 
In this section, we discuss that we can indeed identify the same structure as the U(l) vortices by a simple reduction 
of the U(2) system. 

So far, we have not specified any detail of the interaction, but suppose that one can turn on an interaction which 
breaks the global U(2) symmetry so that only the U(l) symmetry is preserved. Then, in the presence of such an 
interaction, only one Dirac fermion would remain massless. We can realize such a case by simply ignoring the lower 
component ipf {i = 1, - ■ ■ ,n) of the U(2) Dirac fermions. Under this simple procedure, the Hilbert space H{"> of the 
U(2) Dirac vortices reduces to the Hilbert space H'^"^ of the U(l) Dirac vortices discussed in the previous paper [Tl| 
(and also in Sec. Ulland Appendix E)) . Let us discuss the cases of n = 2 and 3, separately. 

Consider first the case of n = 2 (see Sec. lIIIB]) . One finds that only the following states survive after the reduction: 

|0), |2io) - ft^^) |0), |2oi) = f^^+i) |0), |3n) = f"^' o'''' | |0). (4.1) 
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This means that the whole Hilbert space H^^'^ shrinks into its subspaces H-'-i'"'], H^ii"] and H^ini. Notice that these 
four states are equivalent to the basis states for H^^'*'', H*^^'^^ and H^^^^^ in the U(l) Dirac vortices. Thus, it is not 
surprising that the exchange matrices in the reduced U(2) Dirac vortices are equivalent to those of the U(l) vortices: 



l[ool _ (2,0) 
2[iol _ ^(2,1) 



(4.2) 
(4.3) 



Similarly, in the case of n = 3, one finds that only the subspaces H^iffoi^ H^iiffi, H^in"! and H^imi survive after 
the reduction of the lower component. These Hilbert subspaces correspond, respectively, to H^'^'^^ H'^'^'^^ H^'^'^'' and 
][j(3,3) ^jj^g U(l) Dirac vortices. The matrices of these two cases are equivalent to each other: 

T^'-"' = Tf^'\ (4.6) 

r,^'""' = Tf^'\ (4.7) 
r/'"^> = rf^'\ (4.8) 

with £ = fc, fc + 1. 

By repeating the same arguments for general n, one can easily find that there will be a correspondence between 
the Hilbert subspaces H^i^ "i and IHI^"'.'^). The matrices in these two subspaces will be equivalent; 

^A^[i...io...o, ^^(n,/)^ (4.9) 

with €=1,-- - ,ri— 1. Here, the U(2) representation A^n . . . l Q • • • 0] contains / Dirac fermions, and the dimension of 



A^[i...io- -o] is equal to / + 1. Thus, the Hilbert space H^"^ of the U(2) Dirac vortices covers H(") of the U(l) Dir ac 
vortices. This is confirmed also for n = 4 vortices in Appendix [BJ 



V. GENUINE NON-ABELIAN STATISTICS 

A. Sectors of exchanges of the identical vortices 

We have constructed the braid group made of exchanges of U(2) Dirac vortices. One may be able to regard it as 
a kind of exchange statistics since all the Dirac vortices to be exchanged are energetically the same. However, this 
is true only if we are allowed to neglect the difference between the occupation and absence of Dirac fermions in each 
vortex. A genuine non-Abelian statistics in a strict sense should appear in exchanges of identical particles, which 
implies that we should distinguish the number of Dirac fermions in each vortex. In the case of vortices with U(l) Dirac 
fermions, the whole braid group representation is non-Abelian, but does not have a genuine non-Abelian statistics 
when restricted to a subspace of exchanges of identical states (i.e., between two occupied states or two un-occupied 
states). What is truly remarkable in the U(2) case is that the braid group which we have found indeed contains a 
genuine non-Abelian statistics as its subgroup where we consider only the exchanges between the identical vortices, 
the vortices with the same numbers of Dirac fermions. Below, we explicitly show them in subspaces with n = i and 
4 vortices. 



1. The doublet Hilbert subspace H^ii"! for three U(2) Dirac vortices {n = 3; fc-th, (fc + l)-th and (fc + 2)-th 
vortices) with single Dirac fcrmion occupations. The minimum genuine non-Abelian statistics appears in the 
doublet Hilbert subspace H^ii"! in n = 3 vortices. The two-dimensional Hilbert subspace H'^imi is spanned by 
the basis |2aiii) and |2siii) defined in Eq. (|3.31l) . The Dirac fermion numbers in each U(2) Dirac vortex are 

one so that they are all identical vortices. Two exchange matrices, i.e., t^'"^' for the exchange of the k-th and 
{k + l)-th vortices, and t^{!"' for the exchange of the {k + l)-th and {k + 2)-th vortices, can be found in the 
last line in Eq. (|335| : 
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They are non-commutative: '''fc 7^ ''"fc ' showing a genuine non-Abehan statistics. 

2. The non-Abehan transformation equivalent to Eq. (|5.10p are embedded in a sector of four vortices. For instance, 
let us consider four vortices three of which contain one Dirac fermion at each vortex but the rest of which contains 
no Dirac fermion, then the exchange of the identical vortices with the one Dirac fermion occupation leads to the 
same result as in Eq. (|5.10|) . To see this more explicitly, consider the doublet subspacc with the basis |2aoiii) 
and |2som) defined in Eq. (|B.13p . We can see that the submatrices made of the fifth and sixth columns and 

rows in t^{_""^' and t^'^^''"^ in Eq. (|B.36p are nothing but the matrices in Eq. (|5.10p . The rest corresponding 

submatrix in t^'^^^"' in Eq. (jB.36[) is consistently zero, because we are considering the exchange of only three 
identical vortices with one Dirac fermion occupations. 

3. The singlet Hilbert space Hiimii for four U(2) Dirac vortices (n = 4; fc-th, (fc + l)-th, (fc + 2)-th and (fc + 3)-th 
vortices) with single Dirac fermion occupations. We have now four identical vortices each of which contains one 
Dirac fermion. The singlet Hilbert space H-*-!""! is a two-dimensional subspace with the basis |1aaiiii) and 
|lssii n) defined in Eq. (|B.8|) . The non-Abelian exchange matrices can be found in Eq. ()B.3ip as 

i[iiii] _ / -1 \ i[iiii] _ ( h. ^\ - ( ^\ (^\\\ 

~ y I J ' k+i ^ \ V3 _i I ' ^fc+2 ~ y 1 y ' ^ ' 

The exchange of the Dirac vortices are non- Abehan; r/'^^V/jY"' ^ r^j^^^;'""' for ^ = fc, fc + 1. 

4. The triplet Hilbert space H^imu for four U(2) Dirac vortices (ri = 4; fc-th, (fc + l)-th, (fc + 2)-th and (fc + 3)-th 
vortices) with single Dirac fermion occupations. The four identical Dirac vortices with single Dirac fermion 
occupations allow for another sector having a genuine non-Abelian statistics. The triplet Hilbert space H^i""! 
is three dimensional with the basis |3asiiii), |3saiiii) and |3ssiiii) in Eq. (jB.19p . The exchange matrices 
can be found in Eq. ()B.42[) as 

1 \ 

0-10, (5.12) 

1 y 

The exchange of these Dirac vortices are non-Abelian; r^^j""' 7^ r^^j^^'r^^'""' for £ = fc, fc -I- 1. 

From the above examples, we find that the U(2) Dirac vortices with one Dirac fermion occupation give non-Abelian 
quantum statistics of identical particles. 



[1111] 




•i[iiii] 
'fc+i 




\ V2 V2 



[1111] 



fc+2 



B. Entanglement in the Hilbert space 

If one reminds of the case of Majorana vortices showing the genuine non-Abelian statistics, one may critically 
wonder why we have found the genuine non-Abelian statistics even though Dirac fermions are defined locally at each 
vortex. As is well known, for vortices with Majorana fermions, Dirac fermions have to be defined non-locally by using 
two Majorana fermions localized at spatially separated vortices, and, thus far, it has been commonly thought that 
emergence of non-Abelian statistics is attributed to such a non-locality of Dirac fermions. In contrast, in our case, 
Dirac fermions are introduced from the beginning, and are of course defined locally at each vortex. Then, it is natural 
to ask why we get non-Abelian statistics from local Dirac fermions, or to raise question if the non-locality is really 
essential for the emergence of non-Abelian statistics. 

Our answer is that the non-locality is indeed essential but it docs not come from the definition of Dirac fermions 
but rather from the definition of the basis of the Hilbert space. Namely, the use of entangled states, which was 
absent in the previous U(l) Dirac case, is essential in the U(2) case. For instance, in the first example of the previous 
subsection, the basis |2AjLii) and |2siii) in Eq. p. 311) are entangled because of (anti-)symmetric combination of 
indices. Such entangled states are necessary because the basis of the Hilbert space have to belong to the irreducible 
representations of U(2) symmetry. Therefore, genuine non-Abelian statistics docs not appear in U(l) Dirac vortices 
but it first appears in U(2) Dirac vortices because of non-Abelian U(2) symmetry group acting on the doublet Dirac 
fermions on the vortices. 

From these considerations, we conclude that some kind of non-locality is necessary to obtain the genuine non- 
Abelian statistics: It is the non-local definition of the Dirac fermions in the case of Majorana vortices, while it is the 
entanglement in the Hilbert space in the case of Dirac vortices. 
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VI. SUMMARY 

We have considered the simplest system of non-Abehan Dirac vortices, namely the system of vortices in which U(2) 
doublet zero-energy Dirac fermions are trapped. We have constructed the non-Abelian representation of the braid 
group for the exchange of vortices with U(2) Dirac fermions. This was confirmed both at the operator level and in 
matrix representations of the exchange operations. In particular, the whole Hilbert spaces for n = 2, 3, 4 vortices are 
decomposed into subspaces according to representations of U(2), and the matrix forms of the exchange operations 
were obtained in these subspaces. We have found they have off-diagonal elements in some subspaces. This analysis 
is an extension of the previous one for the U(l) (one component) Dirac vortices [ll|. By using a simple reduction of 
the U(2) results, we have indeed identified the same matrices as in the U(l) Dirac vortices. This way of identifying 
the U(l) structure is rather different from the one in the non-Abelian Majorana vortices where one can extract the 
U(l) part as a tensor product. It is not trivial at this point if a similar identification is possible in the Dirac vortices. 
Finally, we have found that the whole braid group contains a subgroup of a genuine non-Abelian statistics for the 
exchange of the vortices with one Dirac fcrmion. The reason why non-Abelian statistics appears only from locally 
defined Dirac fermions is that the basis of the Hilbert space are entangled because of the representation of U(2) acting 
on the Dirac fermions. In this sense, the 'spatial' non-locality of Dirac fermions is not needed to have non-Abelian 
statistics but rather the non-locality, namely entanglement, in representations is needed. 





Majorana vortices 


Dirac vortices 






ik = fk 




# of trapped fermions 


1 


3, 2N-1 


1 


2 


and symmetry 


non 


S0(3), SO(2iV-l) 


U(l) 


U(2) 


examples 


p-wave SC, etc 


CFL phase of color SC 


A-phase of '^He, etc 


CFL phase of color SC, etc 


dimension of 


2?i/2 


j-23yi/2 ^2'^N-lyi/2 


2" 


(2^)" 


n- vortex system 


for even n 


for even n 






representation 


non-Abelian 


non-Abelian 


non-Abelian 


non-Abelian 


of braid group 


(Ivanov matrices) 


(Ivanov X Coxeter) 






exchange 


non-Abelian 


non-Abelian 


Abelian 


non-Abelian 


statistics 


(Ivanov matrices) 


(Ivanov X Coxeter) 






references 


m 


S0(3) [9|, S0(2iV- 1) [10] 


m 


this paper 



TABLE I: Summary of the results on non-Abelian representations of the braid group and non-Abelian statistics of Majo- 
rana/Dirac vortices with Abelian/non- Abelian symmetry. Note that the braid group acts on all particles while the non-Abelian 
statistics concerns only on identical particles. Multiple zero-energy fermions are in the vector representations of the groups 
shown in the table. In the Majorana vortices, the exchange matrices for a single fermion are called Ivanov matrices, and those 
for multiple fermions are decomposed into the Ivanov matrices and generators of the Coxeter group. 



In Table I, we have summarized the results of Majorana/ Dirac vortices with Abclian/non- Abelian symmetry. All 
of them exhibit non-Abelian representations of the braid group. The right most column corresponds to the present 
result. The exchange of Majorana fermions automatically gives non-Abelian statistics, while the exchange of Dirac 
fermions contains non-Abelian statistics when one exchanges identical particles, that is possible at least for U(2). 

So far, we know only high-energy-physics examples of vortices with multiple Dirac fermions , i.e., non-Abelian 
vortices in dense QCD and supersymmetric QCD. As future studies, examples in condensed matter physics and 
experimental observation of the non-Abelian statistics in laboratory will be important. Applications to quantum 
computing should also be explored. 

Appendix A: Exchange matrices for U(l) Dirac vortices 

In this Appendix, we present the matrix representations of the exchange operators t^'s (i = 1, • • • , n — 1) for the 
U(l) Dirac vortices which was first obtained in Ref. [ll[ (below we omit the superscript "s" for notational simplicity). 
The Hilbert space of n vortices, 11^"^ is decomposed into subspaces M^"'-''^ which are specified by the total fermion 
number /. 
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In the case of n = 2 (fc-th and {k + l)-th vortices), we have only one exchange operation T^. The Hilbcrt space 
H(2) is decomposed into a direct sum of three subspaces: H(2^°) = {|00)}, H^^^^) = {|10),|01)} and H(2.2) ^ {|11)}. 
Thus, 11(2) ^ ]g(2,o) Q ]Hl(24) q ]H[(2^2)_ xhen, we obtain the matrix representations of ffc as 

(A.l) 
(A.2) 
(A.3) 

In the case of n = 3 (fc-th, (fc + l)-th and (fc+2)-th vortices), we have two operations Tk and Tk+i. The Hilbcrt space 
H(3) is decomposed into a direct sum of tf^.o) = {|000)}, rf^-i) = {|100), |010), |001)}, H(3^2) ^ {|110), |011), |101)} 
and H(3^3) =^ Namely, H^^) = H(3^°) ® H^^.i) e H(3.2) H^^.a). We thus have the matrix representations of h 

and ffc+i as 



_(2.0) 


= 1, 




_(2.1) 




-1 

1 


_(2.2) 

^k 


= 1. 





J3,0) 



J3,l) 



J3,2) 



J3,3) 



J3,0) 




^(3,1) 
'fc+1 



^(3,2) 
'fc+1 




(A.4) 
(A.5) 

(A.6) 
(A.7) 



In the case of n = 4 (fc-th, (fc + l)-th, {k + 2)-th and (fc + 3)-th vortices), we have three operations Tk, 
Tk+i and Tk+2- The Hilbcrt space H'^'*' is decomposed into a direct sum of five sectors: H^^^o) ^ {|0000)}, 
H(4.i) = {|1000), 10100), 10010), 10001)}, H(4^2) ^ {|1100), |1010), |1001), |0110), |0101), |0011)}, H(4^3) = {|lllO), 
|1101), |1011), 10111)} and H(4.4) = {|1111)}. Namely, U'-'^^ = H'^^") ® Mt^^) ^(4:2) ^ jj(4,3) ^ ]fj(4,4). The ma- 
trix representations of Tk, t^+i and are 



J4,0) 



J4,l) 



(4,2) 



J4,3) 



J4,4) 



_(4,0) _ (4,0) 
'fc+1 ~ 'fc+2 



1, 





(0 




1 


o\ 








(1 














f 


1 





^ 
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(A.8) 



(A.9) 



(A.IO) 



(A.ll) 



(A.12) 



Appendix B: The case of n = 4 for the U(2) Dirac vortices 



We consider the system of n = 4 vortices, which is an ensemble of fc-th, {k + l)-th, (fc + 2)-th, and (fc + 3)-th vortices. 
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The Hilbert space is constructed straightforwardly from the n = 2 and n = 3 cases. The representations of the 
U(2) symmetry are decomposed as 

(lo + 2i + l2)vortcxfc <8) (lo + 2i + l2)vortexfe+l ® (lo + 2i + l2)vortexfc+2 <8l (lo + 2i + l2)vortexfe+3 
= loOOO + l0020 + l0002 + I2OOO + lo200 

+ looii + liloo + liolo + liool + loiio + loioi 

+ I0022 + I2020 + I2OO2 + I022O + I0202 + I22OO 

+ I1120 + I1102 + I20II + I02II + I102I + I1012 + I012I + I0II2 + I2IIO + I2IOI + I12IO + I12OI 

+ I2022 + I0222 + I2220 + I2202 

+ IaA 11 jLl + IsS 11 11 

+ I1122 + I22II + I212I + I2II2 + I122I + I1212 

+ I2222 

+ 2iooo + 2oioo + 2ooio + 2oooi 

+ 2oo21 + 2o012 + 22010 + 22001 + 2o210 + 2o201 + 2i020 + 2i002 + 2oi20 + 2oi02 + 22100 + 2i200 

+ 2aiiio + 2aiioi + 2aioii + 2s ion + 2aoiii + 2soiii + 2siiio + 2siioi 

+ 22021 + 22012 + 2o221 + 2o212 + 22210 + 22201 + 2i022 + 2oi22 + 22120 + 22102 + 2i220 + 2i202 

+ 2ai121 + 2ai112 + 2a2111 + 2s2111 + 2a1211 + 2s 1211 + 2sil21 + 2sill2 

+ 22221 + 22212 + 22122 + 2i222 

+ 3iioo + 3ioio + 3iooi + 3oiio + 3oioi + 3ooii 

+ 32011 + 3o211 + 3i021 + 3ioi2 + 3oi21 + 3oil2 + 32110 + 32101 + 3i210 + 3i201 + 3ii20 + 3ii02 

+ 3as u u + 3sAn 11 + 3ss 11 11 

+ 32211 + 32121 + 32112 + 3i221 + 3i212 + 3ii22 

+ 4s 1110 + 4s 1101 + 4s ion + 4s 0111 

+ 4s 2111 + 4s 1211 + 4s 1121 + 4s 1112 

+ Ssiin, (B.l) 

where the subscript nfenfc_|_ini;+2'T'fc+3 [n-k, ^fc+ij '^fc+2j i^k+s = 0, 1, 2) denotes the number of the Dirac fermions, rik, 
Uk+i nk+2 and Uk+s, at the fc-th, (fc + l)-th, (fc + 2)-th and (fc + 3)-th vortices, respectively. One obtains the basis 
states of the Hilbert space by applying -0"^ {£ — k, k + 1, fc + 2, fc + 3 and a = 1, 2) successively to the Fock vacuum 

|0) defined by i>f\0) = for all £ and a = 1, 2. 
For singlet, there are nine Hilbert spaces, H^poofi, ]HI^[2oooi^ H^inoo]^ M^i^aooj^ H^pnoi, H^[222oi^ H^[iiiii, H^P^ni and 
[22221. They are defined as M^m™"! = {|loooo)} with 

lloooo) = |0), (B.2) 

Hi[2oooi ^ {|loo2o), II0002), II2000), II0200)} with 

|l0020> = ^fc+2^fc+2|0>' 

|looo2> = ^'fc+aV^fc+alO)' 
|l2ooo> ^ i'l^^VlO), 

|lo2oo> ^ ^'ili^'f+ilO), (B.3) 
Hi[iiooi ^ lllooii), lliloo), lliolo), lliool), lloiio), lloioi)} with 



loon) 


— ^\'f^k+2'yk+3 




liioo) 


— ^^^k 'Pk+l 


^r^i;i)io), 


lioio) 




^rv3r+2)io), 


liooi) 
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(B.4) 



•■[2200] — 



{|1 



0022, 



2020/ 



L0220/ 



L0202/ 

II0022) 
II2020) 
II2OO2) 
|lo22o) 
II0202) 
II2200) 



[2110] 



L1120, 



1102/ 



L102I/ 



1 

1 
1 



Oil), |l0211 
|lll2o) 
II1IO2) 
II2OI1) 
|l021l) 
|ll02l) 
II1OI2) 
|l012l) 
II0II2) 
II2II0) 
II2IO1) 
|ll21o) 
|ll20l) 

IHll[2220] ^ {II2022), |1 0222), II2220), II2202)} with 



II2200)} with 

= V'r^f^r+2^'i2io), 

= V'fc+i^'fc+iV'fcU^fcUlo)' 

= V-r+i^'feliV-i+a^fc+slO)' 

= V-r^f^r+i^'lilo), 

111012)7 |loi2l), II0II2), II2II0) 



(B.5) 



L2101/ 



L1210/ 



|li20i)} with 



— (V'r^fcii 



V'fe'v^fcVi)^fe;2^fc;2io), 



It J2t 



fc+3 



It J2t 



1 

1 

71 
1 

71 
1 

71 
1 

71 

^^'r^r(v^r+i^^i2-v3r+i^;u)io)' 



(V^fe+iV'feU - V^r+i^i+3)^fc+2^fcl2lo)' 



(V^fe+iV'feU - V^r+l^il2)^fcl3V'fc+3|0): 



1 

7s 
1 

V2 
1 

72 



>'ll^^l3)|0): 



II2022) 
II0222) 
II2220) 
II2202) 



= V^fe+l^fc+l'/'fc+2V'fc+2^fc+3'^r+3|0)' 



Mini] — 



{|1aa 



11 11/ 



LSS 11 11 



)} with 



|1aaj 



111/ = 2 



^rv'fc+i)(v^fc+2V'fc+3 



(B.6) 



(B.7) 



(B.8) 



18 



LSS 11 11/ 



1 

71 



'ltJ2t 



;2tjit 



'It J2t 



?2t Jlf 



?2tJ2t ,71t 



'It 



with a notation A (S) for an antisymmetric (symmetric) combination in first and second pair of indices with 
the underhne in the fc-th and (fc + l)-th vortices and in the (fc + 2)-th and (fc + 3)-th vortices. M-'-i^aiii = 
{II1122), II2211), II2121), II2112), |li22i), II1212)} with 



II1122) 
II2211) 
II2121) 
II2112) 

|ll22l) 
II12I2) 



1 



V^f '/'fe+l)V'fcV2'/'fe+2^'fc+3'/'fc+3|0>' 



It J,2t Jit J,2t 



1 



^V^fe^V^f (^fc+lV'fcU - ^fcllV'fcl3)^il2'/'fcl2|0>' 
^V'fc^^f ('/'fc+lV'fc+2 - ^fc+lV'fc+2)^fc+3V^fc+3|0)' 

1 

1 



^/2 



V-f ^fe+3)^fc+lV'r+l^fc+2'^fc+2|0)' 
(V'^^fcU - V'fel2)^fc+lV'fe+i^fc+3V^fc+3|0)' 



and ]HI^P222i = {II2222)} with 



II2222) = i^lHlHl\ii^V+i'^l+2iV+2'^k+ii^'^'^ 



fc+3 1 



(B.9) 



(B.IO) 



For doublet, there are six Hilbert spaces, H^iioooi, H^pioo]^ M^imoi, M^i^^mi, H^pmi and H^i^^^i). They are defined as 
e2[ioooi = {|2iooo>, |2oioo>, |2ooio>, |2oooi>} with 



'1000/ 



'0010/ 



'0001/ 



^ivl 

v-^fc+1 

Vk+2 

''-^fc+3 
''^fc+3 



|0>, 
^ |0), 
^ |0), 

1 10), 



(B.U) 



'^[2100] = 



= {|2oo2i), |2 



0012/, P2010/ 



22001), |2o2io) 

|2o02l) 

^0012) 



,|2 



0201/, [■^1020, 



Vk+2 Vk+2 



If 



, |2l002), |2oi2o), |2oi02) 

|0), 



|22ioo), |2i20o)} with 



'f^k+3 
l-^k+3 



V'fc+2 ] Jit J2t |Q^ 
72i V^fc+3''^fc+3lW, 
V'fcV2 / 



22010) = V^fcVfc 



?2t / ^fc 



If 

fc+2 



2t 
fc+2 



22001) 

2o2io) 

2o20l) 



fc+3 
V^fc+3 



|0), 
|0), 



^Ili^^Vi f 1 10), 
|0), 



V^fc-+l''^fe+l 



''^fe+3 
''^fe+3 
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'[1110] 



{|2 



Ajllio), |2a 1101/ 

|2a nio) 
|2a noi) 
I 2a ion) 



'S 1011/ 



I 2a 



0111/ 



'SOlll/ 



's iiio; 



'S 1101/ 



I 2a 



1011/ 



1 

71^ 



1 

71 




(B.12) 



S 1101 



)} with 



^/2 v^r 



('/'fc+2V'fcl3 



It 



'It 



72t^/lt J,2t 



',lt^/,l1^ 



?2t 



|0), 



t-k+i 
'/^fc+i 

72t Jit 



TjV'fe+iV^fcU^fcU - 75V'fe+i('/'fc+2V4l3 + V^U^fcU) 



'It 



It „7,2t 



/2t Jit 



7gV'fc+i(V4+2V4+3 + V'feUV'fcU) - V l^fc+iV'fe+2V'fc;3 



|0), 



;)3C<i^a2 - ^(cCi + V'r^iii)^^i2 



^('?/3^Tj2t , j;2t?lt ^J2t _ /2j2t 72t ?r 

ye'^'fe ^fc+i ^^fe Vk+i)Vk+2 V 3^fc Vk+iVk 



^k+2 



|0), 



|0>, 



(B.13) 



with a notation A (S) for an antisymmetric (symmetric) combination in a pair of indices with the un- 
derhne in the fc-th and (fc + l)-th vortices, or in the (fc + 2)-th and (fc + 3)-th vortices. 

{|2202l), I22OI2), |2o22l), |2o212), |2221o)7 |2220l)7 |2l022), |2oi22), |2212o)5 |22102), |2l22o)i |2l202)} with 



^[2210] — 



|2202l) 



'2012/ 



'0221/ 



y-^k+a 



|0), 



Vk+2 
Vk+2 



k+l't'k+2^k+2 
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2o212) 
22210) 



'1022/ 



'2120/ 



'1220/ 



'1202/ 



i^k^k'^k+i^k+i 




2t 

fc+l 



V'fc+2^fc+2V'fc+3V'fc+3|0)' 



''^fc+1 
V'fc+l 

Vk+1 



V'i+i^fcliV'feU^fcUi 



V'fc+lV^, 



+ 1 Vfc+lV'fe+s^fc+sl^)' 



(B.14) 



^1111 = {|2aii2i), |2aiii2), |2a21ii), |2s2iii), |2ai21i), |2s 1211), |2s 1121), |2s 1112)} with 

?2tJlt 



It J2t 



'A 1112/ 



|2a 



2111/ 



PS 2111/ 

|2a 1211) 
1 2s 1211) 

|2s 1121) 

|2s 1112 



1 

1 



2tJ,it ^ 

k+l) 



Vk+2 



vk+\ 

v-'fe+l 



fe+2 VA;+3 
If 



fc+3 

V'f+2V'fcU)|0)' 



^3 V^fc+l ''^fc+2 ''^fc+3 
_.?,2t ^J,lt J,2t 
76 



2t Jit 



rP,Vk+l\Vk+2Vk+Z + '/^A;+2'/^fc+3/ V 3 ^fc+1 V^fc+2 ''^^ 



If 



2t 1 V'fc+lV'fe+l(V'fe+2^'fc+3 - V'r+2^fel3)|0)' 



'It J,2t 



2t Jit 



?it 



7gV'fe'^(V'fe+2^fe+3 + V'feV2V^fcl3) 



■j^-^k iA+2^k+3 + '^k+2Vk+3) 



/,2t ^It 



2 Tit 72t 72t 



3 "^fc 
It 



^fe+2V'fe+3 



|0>, 



|0), 



J_fJ,itJ2T , J^2TTiT ^J2T _ /2j;2Tj;2T J 



rit 



^fc+lffc+3 



) ^ 



73'^r'/'fe+l'/^fc+2 



2tJ,lt 



?lt 



V3 



(B.15) 



with a notation A (S) for an antisymmetric (symmetric) combination in a pair of indices with the underUne in the 



fc-th and [k + l)-th vortices, or in the (fc + 2)-th and (fc + 3)-th vortices. 



^[2221] 



{|22 



2)} 



21 



with 



|2222l) 



'2212/ 



'2122/ 



'/'fe^V'f V'fc+l^fe+l^fc+2^fe+2 



|0), 



If 



If 



) ^^U^^UV^iUV^^Ulo), 
"^fe+i / 

It J2t Jit J,2t Jit J2t 



2t V'fe+i^fc+iV'fe+2V'fc+2V'fc+3V'fc+3|0)- 



(B.16) 



For triplet, there are four Hilbert spaces, H^inoQi, M^i^noi, H^iimi and H^i^^ni. 
{I^iioo)) l^ioio)) |3iooi)i |3oiio)i |3oioi)7 l^ooii)} with 



They are defined as H^iiiod = 



3[2110] = 



= {|32oii), |3( 



>1100/ 



'rM'€+^ + ^lH'lU |0), 



V2 

'^k^'4'k+2 

^M'^lU+^'kH'Z2) 
"2tJ2t 



|0>, 



^kt 



k rk+2 



>1001/ 



7f(V'r^fcl3 + V'r^ll3) I io>, 



V2 



>0101/ 



JOOll/ 



71 (^fc+lV'fe+2 + ^fc+lV'fe+2) 
'yk+l'fk+2 

'^fc+lV'fc+s 
75(V'fc+iV^fc+3 + V'fe+iV'fe+3) 
V'fe+iV'fc+3 

'^fc+2^fc+3 



'yk+2Vk+3 



|o), 



(B.17) 



0211/, 101021/, 101012/, |'>0121/, |«J0112 



i,|3 



2110/ 



>2101 



i,|3 



1210/, 101201/, |Oll20 



|3ii02)} with 



V^fc+l''^fe+l 



\k+2 

72^^1+2^1+3+ ^V+2i'k+3) 



|0>, 



V'fe+2'/'fc+3 

Vk+2Vk+3 



71 (V'fe+2'/'fc+3 + '^V+2^V+3) 



|0>, 



''^fe+2'f^fc+3 



>1021/ 



'f^fc rk+3 



Vk+2 Vk+2 



^k+3i^k+3\^)^ 



|0>, 



V'f V'fe+2 



22 



^k+l^k+S 



|3oil2) 



Vk+lVk+2 



|32 



101/ 



>1210/ 



>1201/ 



>1120/ 



Vk+lVk+2 
"2t , J2t 



^2{^ll,rkl2 + i^t;+i^ll2) 1 |0), 



V'fc+lV'fc+2 



2t 



7|('/^fe+lV'fe+3 + ^fe+l^fc+3) 
V^/c+lV-'fc+3 



|0), 



^fV'fc+2 



'Pk+irk+l 



|0), 



k+3 



|3iio2> ^ 



%/2 



'[""1 = {|3asiiii), |3saii u), |3ssii n)} with 



Vk rk+i 



|3asj 



V'fe+2'/^fe+3 



It 



'It ,72t 



-^i^k+2^k+3 



+ r,l2iil3) lo), 



2t j;2t 



|3SA 



11 11/ 



'/^fe+2'/^fc+3 



)|0>, 



(B.18) 



|3ssiiii) = 



/ iV'fc^V'fc+iCV'fcUV'fcla +V'fc+2^fc+3) 

^li^U^^ 7/,2t j;2t 



+ ^k^k+l)^k+2^k+3 



, -2t -2t "It ''f^^ 
:y5'Pfe V^fc+lVfe+2'/^fc+3 ~ TfV'fc ^fc+lV'fe+2V'fc+3 



It 



|0>, 



(B.19) 



\ ifW;iU2t +J,2tj;it ^J,2t ?2t _lJ;2tT2t /Tit 72t ?2t Tit ^ 

\ 2V'f^fe V^fc+l V^fc Vk+i)Vk+2Vk+'3 2'f^k rk+l\Vk+2rk+3 ^ ^k+2l'k+3) 

with a notation A (S) for an antisymmetric (symmetric) combination in first or second pair of indices with 
the underhne in the fc-th and {k + l)-th vortices, and in the (fc + 2)-th and (fc + 3)-th vortices. H^i^^i^i = 

{|3221l), |3212l), 132112), |3l22l), |3l212), |3ll22)} with 



522ii) = 4>immii^iii 



V^fc+2''^fc+3 

'Pk+2'Pk+3 

V^fc+l''^fc+3 
J_(^lt J2t ^J2t 



v^^;3+v^^;iV4.i3) I io)> 



^yrk+i'rk+3 



'Pk+l'Pk+3 



23 



|3l22l) 



SII22; 



75 



If 



Vk+lVk+2 



2t 
fc+1 



72 



■^k^k+^ 
It J2t 

V'fc^V'fc+2 



rk Vk+l 



k+lVk+2> 



|0>, 



V'fc+l'/'fc+lV'fc+2V'fe+2 



|0), 



fe+1 "fk+i ''"fc+a^fe+s 1^)' 



/c+2V^fe+2'/^A:+3^fe+3|0)- 



(B.20) 



For quartet, there are two Hilbert spaces, 
{|4siiio), |4siioi>, |4sioii>, |4soiii>} with 



[1110] and 



1[2111] ^ 



They are defined as 



1[1110] 



*S1110/ = 



|4s 



1101/ 



|4s 



1011/ 



isoiii; 



-^'^k^k+l'^k+2 



i\Hl+ii^k+2 

/ltj,2t 



+ i'Vi'lli)i'll2 



V'f^fc+l)'/'fe+3 



|0), 



73 



fc+3 

^k'4'k+2^k+3 



^V^fc V^fc+2'f^fc+3 ^ \/3^fc y^k+2f^k+3 
J,2t , T2t Tit N 



f '/'fc+2V'fe+3) 



|0), 



V'i+lV^fc+2V^fc+3 



-'/'fcl2V4l3) 



^V4li^fc+2^fc+3 + -ji'i'k+MiU'^'V+s 

^'4'V+Mk+2'i'k+3 + V'fc+2^fe+3) + 7fV'fcllV'fc+2^fe+3 



|0>, 



(B.21) 



and 



1[2111] 



{|4s2iii), |4si2ii), |4sii2i), |4siii2)} with 



*s2iii; 



1 4s 1211) — 



( 



V'fe+l'/'fe+2V'fc+3 



^V^fc+l'/^fe+2V^fc+3 ^ ^'f^fc+lV'/'fc+2rfc+3 ^ Vk+2Vk+3l 



^Vk+\\Vk+2Vk+3 



+ 'Pfc+2V^fc+3/ ^ ^V-fc+l'/^fc+2V^fc+3 



V'fctl^fc+2^fc+3 



|0), 



75V'fe^(^fe+2V'fe+3 +^fe+2^fcl3) 



ygV-'fe l,'/^fe+2V'fe+3 



73 



24 



*S 1121/ 



is 1112/ 



V'fe'V'fcVl)'/'fe+3 
^fc^V'fe+lV'fc+3 



72tJ,2t 



fe+2 
If 



(B.22) 



with a notation A (S) for an antisymmetric (symmetric) combination in a pair of indices with the underline in the 
fc-th and (fc + l)-th vortices, or in the (fc + 2)-th and (fc + 3)-th vortices. 

For quintet, there is one Hilbert space H^[""i which is defined as H^[""i = {|5ssiiii)} with 



|5ss 



11 11/ 



^fc^V'feiiV^fci2V'r 



V6 



k-\-l T^k-\-2 



\ 



k+2Vk+?, 

"2tJ,2t 



\{'iki'V+l + V'fe^^fcil)V'feV2V'fcl3 + 5'/'?V'fe+l(^fci2V'fc+3 + i'k+2^k+3) 



\ 



J 



(B.23) 

with a notation S for a symmetric combination in first and second pair of indices with the underfine in the k-th and 
(fc + l)-th vortices, and in the (fc + 2)-th and (fc + 3)-th vortices. Therefore, the Hilbert space is totally given as a 
direct sum, 



e{«=4} ^ 



jji[oooo] ; 

][2[ioooi ; 

][3[iiooi ; 

l4[iiioi ; 

It5[1111] 



l[i[2oooi ; 

i[2[2iooi : 

jSpiio] ] 

|[4[2111] 



[[i[iioo] ; 

[[2|iiio] 0; 
[l3[iiii] fR ■ 



[[1[2200] ]g[l[21101 H"''P^^°1 ' 
[[2[2210] ]fj2[2111] ]g[2[2221] 
[T3[2211] 



■-[11111 



•■[2211] ]g[l[2222] 



(B.24) 



With these basis states in the Hilbert spaces, the operators fk, fk+i and are expressed as matrices. For singlet, 
the matrices are 



■■[0000] 



■■-[0000] 

'k+1 



-•■[0000] 

^k+2 



(B.25) 



for H-'- [00001 , 
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(B.26) 



for H^poooi 
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'[1111]. Interestingly, we find again the non-Abehan matrices in several Hilbert subspaces. For example, the ma- 



trices T, 
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■■(2000] 
'fc+1 



and Tj,^™"' in the Hilbert subspace 



2000] are non-commutative; t. 



42000] -■■[2000] 
'''f.+i 



[2000] ■■-[2000] 
+1 



for I = fc -I- 1. Therefore, the exchange of the ^-th and [I -I- l)-th vortices (£ = fc, fc -I- 1) induces the non-Abelian 
representation of the braid group in IHI-'"P°°°i. Similarly, the non-Abelian representation of the braid group is realized 



in the following Hilbert subspaces; H [n""! , 



■-[2200] 



■■[2110] 



■[2220] 



Mini] and H^[22iii for singlet. 



*[1000] 



^[2100] 



e2(iiio]^ ]g[2(22io]^ IHI^[2iiii and H^[222i] for doublet, IHI^[ii™i, IHI-^[2iioi, H^'iniu and ]HI-*[22iii for triplet, ]HI*[iii'ii and ]HI^[2iiii 
for quartet. 
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As we have discussed in the text, the U(l) Dirac fcrmions arc embedded in the U(2) Dirac fermions. When tjj'^ 
{£ = fc, fc + 1, fc + 2, fc + 3) are set to zero, the Hilbert subspaces H^poooi, H^i"™]^ H^mooi, H^mioi and H^miii in 
U(2) Dirac vortices coincide with the Hilbert subspaces H^^'^), H^^.i), H(4.2)^ ^^^,3) ^^^^ ]fj(4,4) u(l) Dirac vortices, 
respectively. The matrices between the two are equivalent; 

^l[oooo] ^ ^(4,0) ^g^^^^ 

(B.48) 
(B.49) 
(B.50) 
(B.51) 

with i = k, k + 1, k + 2. 



2[iooo] _ ^(4,1) 

— 

3[iioo] _ (4,2) 

4[iiio] _ (4,3) 

5[iiii] _ (4,4) 

If — ' f 



Appendix C: Subspaces with (rfe)'^ = 1 in U(2) Dirac vortices 

We recall that, for both cases with U(l) and U(2) Dirac fcrmions, four-time exchange of vortices is equivalent to the 
identity; (Tfc)'* = 1. The same relation holds at the operator level: (f|)'* = 1 for U(l) Dirac vortices and (ffe)'* = 1 for 
U(2) Dirac vortices. However, a matrix Tk representing f| or fk happens to satisfy a stronger relation, (rfe)" = 1, in 
some Hilbert subspaces. There, two-time exchange of vortices is equivalent to identity. In this Appendix, we explain 
how such a relation can be satisfied, and check if the representation of the braid group is still non-Abelian. 

First, we consider the case of U(l) Dirac fermions. From Eq. (|2.9p . we find that (f^)^ is expressed in terms of the 
number operator of Dirac fermions in the i-th {£ = k, fc + 1) vortices, If we define = 0, 1 as an expectation 

value of the number operator ijf^ we find that the matrix representing the operator f| yields 

(r,^)2 = (l-2iV,^)(l-2A^|+i). (C.l) 

The right-hand-side reduces to 1 only when (Affe,A^|+i) = (0,0) or (1,1). Therefore, we conclude that the relation 
(r|)^ = 1 holds only when all the vortices are empty or fully occupied: (A^f , • • • , A^^) = (0, • • • , 0) or (1, • • • , 1). In 
both cases, t| = 1 for any £ = I,-- - ,n — 1, and hence it gives just trivial representation of the braid group. In 

fact, as shown in Appendix |Al the matrices for the empty and fully-occupied states are r^^"'*^^ = r^^"'"' = 1 with 
£ = 1,. . . ,n- 1 for n = 2, 3, 4. 

Second, let us consider the case of U(2) Dirac fcrmions. In this case, there are non-Abelian matrices satisfying 
(r^)2 = 1 for any £ = 1, • • • , n — 1. From Eq. p.4p . we find again that (Tfc)^ is expressed in terms of the number 
operator of Dirac fermions in the ^-th (£ = fc, fc + 1) vortices, ^V"" with a = 1, 2. Note that the indices a = 1, 2 of 
the pseudo-spin are introduced. Then, defining Ng = 0, 1 as an expectation value of the number operator ip'^^ipg, we 
find that the matrix representing the operator yields 

(Tfc)2 = (1 - 2Nl)il - 2Nl){l - 2Nl+,){l - 2Nl+,). (C.2) 

The relation (Tfe)^ = 1 is fulfilled by the following combinations, 

{NlNlNl+^.Nl^^) = (0,0,0,0), (1,1,1,1), (1,1,0,0), (0,0,1,1), 

(1,0,1,0), (1,0,0,1), (0,1,1,0), (0,1,0,1). (C.3) 

When we define Ni — + Nf , the above combinations are further rewritten as 

{Nk,Nk+i) - (0,0), (2,2), (2,0), (0,2), (1,1). (C.4) 

Therefore, we conclude that the relation (t/;)^ = 1 is satisfied when Nk + Nk+i is an even number. This conclusion 
is consistent with the expectation from the transformation properties of the Dirac fermions under the operation 
(Tfc)2. Under the two successive exchanges of fc-th an (fc -I- l)-th vortices, the Dirac fermion operators ■0^ and i^lj^i 
are multiplied by —1. If a state is composed of an even number of fc-th and (fc -I- l)-th fermions, the minus signs 
cancel and the state is unchanged under (Tfc)^- Therefore, in order for that the condition {tiY = 1 holds for any 
£=!,■■■ , n — 1, a sum of the Dirac fermion number in the every neighboring vortices, Ni + Ng^i, has to be an even 
number. We note that (A^i, • • • , A^„) = (0, • • • ,0) and (2, • • • ,2) corresponding to the empty state and fully-occupied 
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state, respectively, give just trivial representation of the braid group = 1, like the case of U(l) Dirac fermions as 
discussed above. It is also the case for the state with the highest dimension in pseudo-spin representation. However, 
the other combinations of (A^i, • • • , Nn) induce non-Abelian matrices, namely non-Abelian representation of the braid 
group. 

Let us see examples in U(2) Dirac fermions with n = 3 and 4. 

a) For n = 3, the Hilbcrt subspaces, in which Ne + iV^-i-i is an even number for any £ = I, 2, are 

]H[i[ono]^ H^[2™1, H^P2"i, H^[222]^ tfiiiii and H'^iim. (C.5) 
The matrices in each Hilbert subspace are 

^1,000,^ ^l[.oo,^ ^l[..o,^ ^2[,,„ ^4,,,,^ ^^g^ 

Among them, the matrices t^^°°°^ , r^^'^^^' and r^'^^^' are trivial, because they correspond to the empty state, full- 
occupied state and the state with highest dimension in pseudo-spin representation, respectively. The other matrices 
Te^^°°\ and r|''"' are non-Abehan matrices, and hence the Hilbert subspaces H^i^™!, H^i^^oi and H^imi lead to 

the non-Abelian representation of the braid group satisfying (t^)^ = 1 for ^ = 1, 2. 

b) For n = 4, from Appendix [Bl the Hilbert subspaces, in which Ni -\- Ni^i is an even number for any ^ = 1, 2, 3, 
are 

jjipooo]^ lei^i™™!, H^i""™!, 1I^[222«1, H^i""!, ]HI^[2222i^ H^iiiiii and M^^^^^^K (C.7) 
The matrices in each Hilbert subspace are 

"■[0000] ^1[2000] ^1[2200] ^1[2220] .^.^[lllll ^1[2222] ^^[llll] g^^^J .^.^[iiii] (C8) 



T. 



^l[nnnnl If' 

occupied state and the state with highest dimension in pseudo-spin representation, respectively. The other matrices 



Among them, the matrices t« , -f-^P^^^i and r!"'^^^^' are trivial, because they correspond to the empty state, full- 



^ii.uuuj ^ ^^^i22ouj ^ ^-122201 ^ ^^i^^^^i ^ud t/'""' arc non-Abehan matrices, and hence the Hilbert subspaces H^i^""! , H^i^^ooi ^ 
]H[i[222o]^ jjiiiiii] and H^iimi lead to the non-Abelian representation of the braid group satisfying (rf)^ = 1 for £ = 1, 
2,3. 

We note that the condition (rfe)" = 1 with the braid relations (i) TkT£Tk — TgTkTf for j/c — £| = 1 and (ii) T^Tf ~ tit^ 
for \k ~ i\ > 1 leads to the relations, (i') (rfer^)^ = 1 for |A: - £| = 1 and (ii') {rknY = 1 for |A: - £| > 1 for 
k,i= I,--- ,n— 1. The relations (i') and (ii') imply that such matrices are regarded as generators of the symmetric 
group or the Coxeter group of the type An-i- We recall that the Coxeter group was obtained in the S0(3) or more 
generally SO (A'') (iV odd numbers) symmetric Majorana vortices [1, H^. There, the matrices for exchanging vortices 
were tensor product of "the Ivanov matrices" found by Ivanov Q and the generators of the Coxeter group of ^2m-i 
type {2m the number of vortices). In the present case of U(2) Dirac vortices, however, such tensor structure was not 
found. 
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